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Abstract 

The collision of test charged particles in the vicinity 
of an event horizon of a weakly magnetized non-rotating 
black hole with gravitomagnetic charge has been stud- 
ied. The presence of the external magnetic field de- 
creases the innermost stable circular orbits (ISCO) radii 
of charged particles. The opposite mechanism occurs 
when there is nonvanishing gravitomagnetic charge. 
For a collision of charged particle moving at ISCO and 
the neutral particle falling from infinity the maximal 
collision energy can be decreased by gravitomagnetic 
charge in the presence of external asymptotically uni- 
form magnetic field. 
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1 Introduction 

At present there is no any observational evidence 
for the existenc e of gravitomagnetic monopole, i.e. 



19631 ) parameter or 



so-called NUT (INewman et al 
magnetic mass. Therefore study of the motion of 
the test particles and particle acceleration mecha- 
nisms in NUT spacetime may provide new tool for 
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studying new important general rclativistic effects 
which are associated with nondiagonal components 
of the metric t ensor and have no Newtonian ana - 
logues (See, e.g. Nouri-Zonoz and Lvnden-Belll (1997 ): 



iLvnden-Bell and Nouri-Zonoz ( 19981) :lNouri-Zonoz et. aM 



iNouri-Zonoa (120041) : iKagramanqva et al 



(200c 



2010l) : |m orozova and Ahmedov! |2009l) where solu 



tions for electromagnetic waves and interferometry 
in spacetime with NUT parameter have been stud- 
ied.). Kerr-Taub-NUT spacetime with Maxwell and 
dilation field s is recentl y inves tigated by authors of 
" (120081) . 



the paper 
paper s 



Aliev et al 



Morozova et al 



In our preceding 



Abduiabbarov et al.l 



( 20081) we have studied the plasma magnetosphere 
around a rotating, magnetized neutron star and charged 
particle motion around compact objects immersed in 
external magnetic field in the presence of the NUT 
parameter. The Penrose process in the spacetime of ro- 
tating black hole with nonvanis hing gravitomagnetic 



charge has been considered in lAbduiabbarov et al 



( 201l[) . The electromagnetic field of the relativistic star 
with nonvani shing gravitomagnetic charge has been 
considered by I Ahmedov et. all (|2012l) . 

Astrophysical processes which may produce high en- 
ergy radiation near a rotating black hole horizon at- 
tract more attention in recent publicati ons. Such pro - 
cesses connected with the Penrose e ffect ( Penrose] 1969 ) 
have been studied l o ng tim e ago by Piran et. all dl97oTl: 
Piran and Shahaml (jl977al ibh. Recently iBanados et al.1 



(|2009l) (BSW) demonstrated that for an extremely 
rotating black hole such collisions can produce par- 
tic les of high ce n ter-of- mass energy. The results 
of Banados et al.l d2009l) have been commented by 



Berti et al.1 (|2009l ) where authors concluded that as- 
trophysical limitations on the maximal spin, back- 
reaction effects and sensitivity to the initial condi- 
tions impose severe limits on the likelihood of such 
accelerations. The above mentioned discussions forces 
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us to study the particle acceleration in the space- 
time of a black hole with nonvanishing gravitomag- 
netic charge. The acceleration of particles, circular 
geodesies, accretion disk, and high-energy collisions 
in the Janis- Newman-Winicour sp a cetime have been 



consid ered by iPatil and Joshil (|2012l) ; IChowdhurv et al 
( 20121) . 



The aim of this paper is to show that a similar ef- 
fect of particle collision with high center-of-mass energy 
is also possible when a black hole is non-rotating (or 
slowly rotating) provided there exists magnetic field in 
its exterior and in the presence of the gravitomagnctic 
charge. It is very interesting to study the electromag- 
netic fields and particle motion in NUT space with the 
aim to get new tool for studying new important general 
relativistic effects. This observation might be interest- 
ing since there exist both theoretical and experimental 
indications that such a magnetic field must be present 
in the vicinity of black holes. In what follows we assume 
that this field is weak and its energy-momentum docs 
not modify the background black hole geometry. For a 
black hole of mass M this condition hol ds if the strength 
of magnetic field sa tisfies the condition (jPiotrovich et.al 
201lHFrolovll2012l) 



B rnax 



Mq 

G 3 / 2 M q V M 



10 1 



M 



Gauss . (1) 



We call such black holes weakly magnetized. One can 
expect that the condition (TTJ) is satisfied both for stellar 
mass and supermassive black holes. 

The paper is organized as follows. Section II is de- 
voted to study of electromagnetic field and charged par- 
ticle motion in the magnetized black holes with NUT 
parameter, with the main focus on the properties of 
their ISCOs. Particles collisions in the vicinity of a 
weakly magnetized black hole with nonvanishing grav- 
itomagnctic charge have been studied in Section III. 
The concluding remarks and discussions are presented 
in Section IV. 

Throughout the paper, we use a space-like signature 
(— , +, +, +) and a system of units in which G = 1 = c 
(However, for those expressions with an astrophysical 
application we have written the speed of light explic- 
itly). Greek indices are taken to run from to 3 and 
Latin indices from 1 to 3; covariant derivatives are de- 
noted with a semi-colon and partial derivatives with a 
comma. 



2 Charged particle motion around black hole 
with nonvanishing NUT charge 



external static axisymmetric and uniform at the spatial 
infinity magnetic field. The ap propriate spacetime met- 
ric has the following f orm (see iDadhich and Turakulov 
2002HBini et~ai1l2003l ): 



ds 2 = -— dt 2 +4^ I cos Odtdip + ^ (E 2 sin 2 9 - 

Zj Zj Zj 

- 4A; 2 cos 2 9)d<p 2 + — dr 2 + Y,d9 2 , (2) 

where parameters £ and A are defined as 

Z = r 2 + l 2 , A = r 2 -2Mr-l 2 , 

where I is the gravitomagnetic monopole momentum. 
Here we will exploit the existence in this spacetime of a 
timelike Killing vector £?U = dx a /dt and spacelikc one 
£?n — dx a /dip being responsible for stationarity and 
axial symmetry of geometry, such that they satisfy the 
Killing equations 



(3) 



which gives a right to write the solution of vacuum 
Maxwell equations UA^ = for the vector potential 
of the electromagnetic field in the Lorentz gauge in 
the simple form 



»(¥>)■ 



(4) 



The constant C2 = B/2, where gravitational source is 
immersed in the uniform magnetic field B being paral- 
lel to its axis of rotation. The value of the remaining 
constant C\ can be easily calculated from the asymp- 
totic properties of spacetime ((2} at the infinity. Indeed 
in order to find the remaining constant one can use the 
electrical neutrality of the black hole 

AttQ = ^fF aP 9 dS afi =C 1 f^T a m p q t) {Tk)dS 

+^T^T a m0e w (rk)dS = O, (5) 

evaluating the value of the integral through the spher- 
ical surface at the asymptotic infinity. Here the equal- 
ity = — £a;/3 = — r^^y following from the Killing 
equation was used, and element of an arbitrary 2- 
surface dS a P is represented in the form (see, e.g. 
Ahmedov and Rakhmatovll2003l ) 



(6) 



Here we will consider a charged particle motion in the 
vicinity of a black hole of mass M in the presence of an 
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and the following couples 

™ a y/1 + (rfc) 2 ' 
Ha ~ y/1 + (rfc) 2 ' 

k a = -{rk)T a + y/1 + (rfc) V qp %to p ^ 

are established between the triple {k, m, n} of vectors, 
n a is normal to 2-surface, space-like vector m a belongs 
to the given 2-surface and is orthogonal to the four- 
velocity of observer t°, a unit spacelike four- vector k a 
belongs to the surface and is orthogonal to m a , dS 
is invariant element of surface, A denotes the wedge 
product, * is for the dual element, r/ a p 7 s is the pseudo- 
tensorial expression for the Levi-Civita symbol tap-yS- 

Then one can insert to = —(1 — M/r), m 1 = 
(1 — M/r), and asymptotic values for the Christoffel 
symbols T% = M/r 2 and T a 13 = -1(1 - 2M/r)cos9/r 
in the flux expression ((SJ) and get the value of constant 
C\ = 0. Parameter I does not affect on constant C\ 
because the integral J* cos 9 sin 6d6 = vanishes. 

Thus the 4- vector potential of the electromag- 
netic field will take the following form 



A 

An = Bl COS ( 

E 



A 3 = -E 2 sin 2 9 - 2 A/ 2 cos 2 I 



(7) 
(8) 



The orthonormal components of the electromagnetic 
fields measured by zero angular momentum observers 
(ZAMO) with four velocity components 



(r a ) 



ZAMO 



K 2A2cos6> \ . , 



AEsin 2 #' VAEftsin 2 ! 



(r a ) 



1 AE sin 2 9 



ZAMO 



K 



are given by expressions 



Brl 



E r 



1 



M 



B' 
B" 



K 
Bl fA 

e 2 " \Jn 

E 2 + 
Bt&n9 



Br 



-,0,0,0 



sin26> 



(10) 



(11) 



2A; 2 cos( 



2 cos 9 
sin 2 9 



TZ-Y, 2 ) , 



sin 2 6> ,(12) 
(13) 



A-E 1 



M 
r 



4l 2 + E 2 tan 2 



(14) 



where the following notation has been used: 

K = E 2 sin 2 9 - 4A1 2 cos 2 6. 

Astrophysically it is interesting to know the limiting 
cases of expressions (fTTj) - (fTI)) . for example in either 
linear or quadratic approximation 0(a 2 /r 2 , 1 2 /r 2 ) in or- 
der to give physical interpretation of possible physical 
processes near the slowly rotating relativistic compact 
objects, where they take the following form: 



E r 
E" 
B r 
B° 



2Blcos9 



r 

Bl sin 9 

r 

Bcos9 
Bsm9 



1 - 
1 + 
1 + 2 



3M 

r 

2cos6» 
sin 2 9 
I 2 1 



r 2 sin 2 I 



1 



M 
r 



16r 2 sin 2 9 
x [Al 2 - AM 2 + 4(7Z 2 + M 2 ) cos 29 



(15) 
(16) 
(17) 

(18) 



In the limit of flat spacetime, i.e. for M/r — > and 
l 2 /r 2 -> 0, expressions ([T5" ]) - ([l8]) give 



E r 



E" -> 



B r -> B cos ( 



B° -> B sin t 



(19) 



(20) 



As expected, expressions (|T9|) -(|20 |) coincide with the 
solutions for the homogeneous magnetic field in New- 
tonian spacetime. 

The dynamical equation for a charged particle mo- 
tion can be written as 



dv^ 



(21) 



where r is the proper time, is the particle 4- velocity, 
u^Ufi = —1, q and m are its charge and mass, respec- 
tively. F a p — Ap,a — - A a ,f3 is the antisymmetric ten- 
sor of the electromagnetic field, which has the following 
four independent components 



B 

E 2 " 
B 



F ai = =2i[(r-M)- Ar] cos( 



F 02 = —Al sin 9 , 

*13 = -T^8l 2 [(r - M)E - Ar] cos 2 6 + Br sin 2 i 
E 2 



^ 23 = |(4/ 2 A + iE 2 )sin20. 

While considering the charged particle motion around 
black hole immersed in the magnetic field it is easy 
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to use two conserved quantities associated with the 
Killing vectors: the energy £ > and the generalized 
azimuthal angular momentum C € (— oo,+oo): 



£ 



P 



mA dt 



41 cos 



dtp 



Bl cos 6»)(22) 



L -H<t>) 



Here P u 



P, 



A 

= — Ami — cos/ 
-4Z 2 — cos 2 I 



dr 

- (E sin 2 6» 



is the gener alized 4-momentum 



parameter 56 are neglected. In the Fig. [T] the radial 
dependence of the effective potential of the radial mo- 
tion of the charged particle is presented for the different 
values of dimcnsionlcss gravitomagnetic parameter l/M 
and magnetic parameter b = qBM/m. 

From the Fig. [T] one can obtain the behavior of the 
charged particle motion in the presence of both the 
gravitomagnetic charge and magnetic field. In the pres- 
ence of the gravitomagnetic charge the minimum of the 
effective potential shifts towards the observer at the in- 
finity which means that the orbits of the charged parti- 



of the particle. It was first shown bv lZimmerman and Shahiff s may become unstable. The minimum value of the 
( 19891 ) for s pherical symmetr ic case (NUT spacetime) 



and later bv lBini et al 



(2003) for axial symmetric case 
(Kcrr-Taub-NUT spacetime) that the orbits of the test 
particles are confined to a cone with the opening an- 
gle given by cos0 = 2£l/C. It also follows that in 
this case the equ ations of motion on the cone depend 
on I only via I 2 ( Bini et all 12003 : Abduiabbarov et al. 
20081 ). The main point is that the small value for 



the upper limit for gravitomagnetic moment has been 
obtained by comparing theoretical results with ex- 
perimental data as (i) I < 10~ 24 from the gravita- 
tional microlensing ( Rahvar and Habibi 2004) . (ii) I < 
1.5T0 -18 fro m the interferometry experimen ts on ultra- 
cold atoms ( Morozova and Ahmedov 2009| ). (iii) and 
similar limit has been obtained fro m the experiments 



on M ach-Zehnder interferometer (jKagramanova et al 
20081 ). Due to the smallness of the gravitomagnetic 



charge let us consider the motion in the quasi-equatorial 
plane when the motion in 9 direction changes as 
8 = 7r/2 + 68 (t), where (50(f) is the term of first or- 
der in Z, then it is easy to expand the trigonomet- 
ric functions as sine 1 = 1 - 58 2 {t)/2 + 0(S8 4 (t)) and 
cos0 = <50(t) - £>(<50 3 (t)). Neglecting the small terms 
0(S8 2 (t)), one can easily obtain the geodesic equation 
in the following form 



dt 


£ E 


d^ 


m A ' 


dtp 


C 






dr 


( L - 




\ m 



qB 
2m 



U 



where the effective potential is 
C/=|(1 + S X 2 ) , 



and 



X = 



c 



qB 
2m 



(24) 
(25) 
(26) 

(27) 
(28) 



In the expressions ([2"4")) -(f2"6l the terms being pro- 
portional to the second and higher orders of the small 



radius of the stable circular orbits increases. The influ- 
ence of the magnetic field has the opposite effect: the 
presence of the magnetic field decreases the minimum 
value of the circular orbits radius and charged particles 
may come much closer to the central object. 

In order to study innermost stable circular orbits 
(ISCO) we use its first and second derivatives of the U 
and equalize them to zero: 



U' =- 



U" 



2Ar(l + 2Cx) 



S 2 



= 1 



A 
" E 
2 

+ E 



4Ar 2 _ 



= 



(29) 



3A 



- 2 



E 4 

x(l + 2£ X )+ +|(1 + E X 2 ) =0 



Now we have two equations with three unknown quan- 
tities C, r and x- Solving the equation (|29|) we derive 
C in terms r and x 



C 



-2Ar 2 + AE + E 2 



AE 2 X 2 



s 3 x 2 



4Ar 



'X 



(31) 



Substituting this equation into (|30p one can easily ob- 
tain the equation expressing x hi terms of r 



X±=± 



where 




(32) 



A = 8A 2 r 2 - 5A 2 E + 12Ar 2 E - 8AE 2 - 3E 3 , 
K = -2A 2 r 2 + 2A 2 E - 4Ar 2 E + 3AE 2 + E 3 . 

The signs ± correspond to co-rotating and counter- 
rotating particle orbits, respectively. The dependence 
of the x± horn the ISCO radius are shown in Fig. [2j 
The shift of the shape of x+ to the right direction 
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a) r b) r 

Fig. 1 The radial dependence of the effective potential of radial motion of the charged particle for the different values of 
the gravitomagnetic charge (a): l/M = (solid line), l/M = 0.4 (dot-dashed line), l/M — 0.8 (dashed line) and for the 
different values of the dimensionless magnetic parameter b = qBM/m (b): b = 0.2 (solid line), b = 0.4 (dot-dashed line), 
b = 0.8 (dashed line). 



with increasing the gravitomagnetic charge corresponds 
to increasing of ISCO radius in the presence of NUT 
charge. 

From the Fig. [5]one can see that X- function steadily 
increase when particle approach to the black hole hori- 
zon and go to infinity in the case when gravitomagnetic 
charge vanishes. In the case of nonvanishing NUT- 
parameter, the function x~ f aus m the point of singu- 
larity. But by reason of that jumps are located inside 
the horizon, the effects that take place there are not ob- 
servable relatively to detached observer and cannot be 
interpreted as some full physical theory In other hand 
the spacetime of nonrotating black hole allows us to find 
analytic continuations of our theories inside the black 
hole horizon up to a point of singularity by excepting 
it. 

Consider a circular motion of the charged particle 
in the presence of NUT charge. The effective potential 
has minimum value at the radius of such orbit. The 
4-momentum of a test particle at the circular orbit of 
radius r is 



P' 

e (t) 



m l( e (t)+ ve (4,))' 
(E/A^/^^E/A) 1 /^, 



(33) 
(34) 
(35) 



Here v (which can be both positive and negative) is a 
velocity of the particle with respect to a rest frame, and 
7 is the Lorentz gamma factor. From normalization 
condition p 2 = — m 2 one has 7 = (1 — t> 2 ) -1 / 2 . For 
q > the Lorentz force acting on a particle with v > is 
repulsive (i.e. directed outwards the black hole), while 
for v < it is attractive. 



Using relation d<j)/dr = vy/r and ([25]) with 9 = n/2 
one gets 



r 



A 



This relation allows us to find 



= yl + r 2 x 2 and v = 



r\ 



i/l + r 2 x 2 



(36) 



(37) 



Using (f3"T)) one can find the values of the velocity v± 
and 7-factor 7± . Fig. [3] shows the velocity of a particle 
at the ISCO as a function of its radius, while Fig. [4] 
shows the dependence of j± from t isco . 
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1 2 3 4 5 6 7 
r/M 

Fig. 3 Velocity of the particles at r ISCO as a function of 
its radius for different values of the gravitomagnetic charge: 
l/M = (solid line), l/M = 0.4 (dot-dashed line), l/M = 
0.8 (dashed line). 

The Fig. [3] shows the dependence of the velocity of 
the particle at ISCO. Since there are two values of ve- 
locity for the same radius one can interpret them as 




a) r/M b ) r/M 

Fig. 2 x+ ( a ) an d X- 00 as a function of the ISCO radius for different values of the gravitomagnetic charge: l/M = 
(solid line), l/M = 0.4 (dot-dashed line), i/Af = 0.8 (dashed line). 




Fig. 4 7+ (a) and 7_ (b) as a function of the ISCO radius for different values of the gravitomagnetic charge: l/M = 
(solid line), l/M = 0.4 (dot-dashed line), l/M = 0.8 (dashed line). 



two values of velocity which correspond for two oppo- 
site directions of motion of the particles. Since charged 
particle and magnetic field interaction depends on the 
velocity direction there should be two values v± for each 
r isco ■ Furthermore the absence of external magnetic 
field (right border of the plots) one can obtain only one 
solution for the velocity at ISCO. One should mention 
that in the non-relativistic case one can get the Keple- 
rian velocity profile. 

In the next step we will consider the case, when a 
neutral particle freely falling from infinity collides with 
a charged particle revolving at the circular orbit near a 
weakly magnetized black hole with nonvanishing gravit- 
omagnetic charge. One can denote by p the momentum 
of this particle, and by m and q its mass and charge. 
Denote by [i the mass of a freely falling particle, and 
by k its 4-momentum. At the moment of collision the 
4-momentum is 

P = P + k, (38) 



and the corresponding center-of-mass energy E c , m , is 

Ec. m . =m 2 +M 2 -2(p,k). (39) 

Using the equation of particle motion around black 
hole with nonvanishing gravitomagnetic charge (1241) 
(|2l)|) one can obtain the following relation for the center- 
of-mass energy of two colliding particles for the weak 
magnetic field approximation: 

^ * O^/^EE*/' . (40) 

In the Table Q] the dependence of the ISCO radius 
and the center of mass energy of colliding charged par- 
ticles have peen shown. From the results on can con- 
clude that gravitomagnetic charge correction prevents 
the particle from the infinite acceleration. 

From the obtained result one can observe that the 
presence of the gravitomagnetic charge will decrease the 
value of the center of mass energy. The role of the mag- 
netic field in particle accelerating process is to decrease 
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Table 1 The dependence of the center of mass energy and ISCO radii from the magnetic parameter b for the different 
values of the specific gravitomagnetic charge l/M: 



l/M 0.2 0.4 0.6 0.8 1.0 4.52 



2 + 0.58&- 1 2.02 + 0.586" 1 2.08 + 0.596" 1 2.17 + 0.606- 1 2.28 + 0.616" 1 2.41 + 0.62&- 1 5.63 + 0.746" 1 



E c . m ./m 1.7476 1 / 4 1.7456 1 / 4 1.7386 1 / 4 1.725& 1 / 4 1.7086 1 / 4 1.6886 1 / 4 1.129& 1 / 4 



the innermost stable circular orbits radii. As particles 
come closer to the black hole horizon their energy at in- 
finity is going to increase. The role of the gravitomag- 
netic charge in this process is opposite: the presence of 
the gravitomagnetic charge increase the radii of ISCO. 



3 Conclusion 



In this report we have obtained the expressions for the 
energy and angular momentum as well as ISCO of the 
charged particle in the vicinity of the black hole in pres- 
ence of gravitomagnetic charge and exterior magnetic 
field. 

Recently Banados et al. ( 20091 ) underlined that a ro- 
tating black hole can, in principle, accelerate the parti- 
cles fal ling to the cen tral black hole to arbitrary high en- 



ergies. iFrolovl (|2012l ) has shown that the magnetic field 



could play a role of charged particle accelerator near the 
nonrotating black hole. Because of some mechanisms 
such as astrophysical limitations on the maximum spin, 
backreaction effects, upper limit for magnetic field, and 
sensitivity to the initial conditions, there appears to be 
some upper limit for the center of mass energy of the in- 
falling particles. One of the mechanisms offered in this 
paper is appearing due to the gravitomagnetic charge 
correction which prevents the particle from the infinite 
acceleration. 
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